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ABSTRACT 
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^  A  method  for  analyzing  stresses  and  displacements  in  filament- 
wound  shells  of  revolution  is  presented.  This  method  of  analysis  is  based 
on  the  membrane  theory  of  shells  but  incorporates  the  edge  effect  theory 
to  predict  a  more  realistic  behavior  in  the  vicinity  of  polar  openings.  In 
addition,  because  it  includes  the  contribution  of  the  matrix,  the  present 
membrane  theory  is  more  accurate  for  the  analysis  of  shells  under 
nondestructive  operating  loads  than  is  the  so-called  netting  analysis  which 
considers  the  filaments  to  be  the  only  load- carrying  members.  The 
membrane  theory  is  also  modified  to  give  accurate  solutions  in  the  shallow 
region  of  composite  shells  and  indicates  significant  bending  stresses  due 
to  orthotropy  and  the  shape  of  the  meridian.  While  these  stresses  are 
localized  near  the  apex  in  thin  shells,  they  may  extend  well  down  along  the 
meridian  in  thick  shells.  / 
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MEMBRANE  AND  EDGE  EFFECT  SOLUTIONS  OF 
FILAMENT-WOUND  SHELLS  OF  REVOLUTION 


Introduction 


Filament-wound  shells  of  revolution  are  widely  used  in  modern  tech¬ 
nology  as  applied  to  aircraft,  missile  and  spacecraft  structures,  pressure 
vessels,  deep  submergence  vehicles  and  machine  construction.  They  are 
formed  by  winding  resin-coated  high-strength  filaments  onto  a  mandrel  and 
then  heat-treating  the  resulting  shell  structure  to  produce  a  rigid  composite. 

The  filaments  are  deposited  in  pairs  that  are  symmetric  with  respect 
to  a  meridian,  along  lines  defined  by  the  intersection  of  the  surface  of  the 
shell  with  a  plane  making  an  angle  y  with  the  axis  of  revolution  z.  Such  a 
winding  pattern  is  referred  to  as  planar-wound  filamentary  pattern.  This 
technique  of  winding  has  found  wide  acceptance  in  the  industry  because  it  is 
easy  to  achieve  in  practice. 

The  present  report  describes  an  analytical  method  of  stress  analysis 
of  filament-wound  shells  of  revolution.  This  method  of  analysis  is  based  on 
the  membrane  theory  of  shells  but  incorporates  the  edge  effect  theory  to 
predict  more  realistic  stresses  in  the  vicinity  of  the  polar  opening  where 
reinforcement  devices  and  the  shell,  each  of  different  elastic  properties, 
interact.  The  local  stresses  near  the  line  of  distortion  (i.  e. ,  the  supported 
edge  of  the  shell)  are  formulated  on  the  basis  of  Gol'denveizer's  edge  effect 
theory  (Ref.  1)  and  then  subsequently  reformulated  with  the  use  of  Steele's 
class  notes  (Ref.  2)  based  on  the  method  of  asymptotic  integration  developed 
by  Blumenthal  (Ref.  3),  for  the  purpose  of  correlating  the  consistency  of 
the  two  methods.  The  two  methods  yield  identical  results  if  only  the  leading 
term  of  asymptotic  solution  is  used. 

A  computer  code  based  on  the  present  theory  has  been  developed 
A  user's  manual  for  the  program  including  a  study  of  the  behavior  of  the 
filament-wound  spherical  shells  will  be  published  at  a  later  date.  The 
code  will  automatically  plot  the  variation  of  the  stresses  along  the  meridian 
as  well  as  through  the  thickness  of  the  shell  in  the  filament,  matrix,  and 
liner. 
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Region  of  Applicability  of  Membrane  Theory 


The  membrane  theory  gives  an  approximately  correct  picture  of  the 
state  of  stress  and  deformation  of  a  shell  only  at  sufficient  distances  away 
from  the  lines  of  distortion  of  the  state  of  stress  on  the  middle  surface. 

The  lines  of  distortion  may  be:  1)  the  edges  of  the  shell;  2)  lines  along 
which  discontinuities  of  the  components  of  the  external  surface  loads  occur; 

3)  lines  along  which  the  middle  surface  of  the  shell  has  an  abrupt  change  in 
curvature;  4)  lines  along  which  the  rigidity  of  the  shell  or  its  thickness  under¬ 
goes  sudden  changes. 

When  treating  a  shell  which  is  not  locally  spherical  and  isotropic  in  the 
neighborhood  of  the  apex,  the  membrane  solutions  are  not  generally  valid. 
Recent  work  by  C.  R.  Steele  (Ref.  4  )  indicates  that  for  orthotropic  shells 
in  the  shallow  region  the  particular  solutions  to  the  shell  equations  do  not 
approach  the  "membrane"  solutions.  It  is  interesting  to  note  that  the  par¬ 
ticular  solutions  show  that  bending  stresses  will  be  present.  These  bending 
stresses  are  not  due  to  edge  effect  solutions,  but  exist  simply  because  of 
orthotropy  and  shape  of  the  meridian  curve. 


Stress -Strain  Relations  for  an  Anisotropic  Elastic  Plate 


The  stress -strain  relations  for  a  general  anisotropic  material  is  of 
the  form  (Ref.  5) 
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where  S..  are  the  "moduli  of  compliance";  and  C..  are  the  "moduli  of  elas- 
13  U 

ticity"  of  the  anisotropic  plate.  They  are  related  by  the  following  expressions 
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If  the  lamina  is  orthotropic  and  the  principal  elastic  axes  coincide  with 
the  coordinate  axes,  =  S^g  =  0 

and  the  moduli  of  compliance  S. .  can  be  expressed  in  terms  of  the  elastic 
moduli  and  Poisson's  ratio  as  follows 


.0  =  c 

E  ’  12  21 

x 


V 

E 


X 


V 

E 

y 


(3a) 


Membrane  Solution 


The  membrane  theory  formulated  in  this  section  is  based  on  the  assump 
tion  that  both  the  resinous  matrix  and  the  filaments  contribute  to  the  strength 
of  the  composite  shell.  This  present  membrane  theory  is  more  realistic 
than  the  so-called  netting  analysis  widely  used  in  the  industry,  which  is  a 
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limiting  case  of  the  present  theory.  The  netting  analysis,  based  on  the 
assumption  that  the  filaments  are  the  only  load  carrying  components  of  the 
composite  shell,  is  inadequate  for  furnishing  the  stresses  developed  in  the 
shell  under  nondestructive  operating  loads.  Neither  is  it  capable  of 
permitting  determination  of  the  membrane  deformation,  nor  of  permitting 
incorporation  of  refinements  necessary  for  handling  the  shallow  region  of 
the  shell  because  of  orthotropic  construction  (Ref.  4),  nor  of  permitting 
the  formulation  of  edge  effects. 


For  a  shell  of  revolution  subjected  to  loads  having  axial  symmetry,  the 
stress  resultants  and  N^  in  meridional  and  circumferential  directions, 
respectively  (Figure  1),  are 
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For  uniformly  distributed  loading  on  a  shell  without  polar  opening,  Eqs.  (4) 
reduce  to 
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where  the  principal  radii  of  curvature  r^  and  r^  are  given  by  the  formula 
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when  the  meridian  curve  of  the  shell  of  revolution  is  defined  by  a  single 
equation 


r  =  f(z) 


(6) 
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Figure  1,  Stress  Resultants  and  Applied  Loads 
on  Differential  Element 
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If  the  shell  is  described  by  a  pair  of  parametric  equations 


r  =  r(s) 
z  =  z{s) 


(7) 


the  radii  of  curvature  can  then  be  computed  from  the  following  relationships 
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where  the  primes  refer  to  the  derivatives  with  respect  to  the  parameter  s 
which  is  defined  herein  as  the  arclength  of  the  meridian  curve. 


For  the  given  stress  resultants 
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where  0 ,  0  are  associated  with  x,  y,  respectively. 


Since  the  filament  patterns  of  composite  shells  of  revolution  are 
generated  by  the  superposition  of  n/2  different  pairs  of  layers  which  are 
symmetric  with  respect  to  the  meridians,  with  2A.  representing  the  cross- 
sectional  area  of  the  filaments  of  one  distinct  pair  and  d.  the  spacing  of  the 
successive  filaments,  the  ratio 

A. 

1 


d. 

1 


represents  the  thickness  a^  of  an  equivalent  solid  sheet  of  filament  (Figure  2). 
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Typical  Meridian 
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Figure  2.  Element  of  Shell  Illustrating 
Filament  Pattern 


Assuming  that  each  filament  crosses  the  equator  and  all  other  parallel 
circles  of  the  shell  an  equal  number  of  times,  the  spacing  of  the  filaments 
along  a  parallel  circle  corresponding  to  a  particular  value  of  the  position 
angle  (t>  can  be  expressed  as 
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where,  N  is  the  number  of  circuits  of  filaments  around  the  shell. 

In  order  to  present  a  method  of  analysis  of  greater  generality,  no  specific 

filament  pattern  is  assumed  and  thus  the  filament  layers  are  variously  oriented 

with  reference  to  a.. 

1 

Hence 
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where  is  the  angle  that  the  filaments  of  the  i-th  layer  makes  with  the 
parallel  circle. 

The  total  thickness  of  equivalent  solid  sheets  of  filament  corresponding 
to  a  particular  (j)  is 
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Hence,  the  net  matrix  thickness 


'’m  '  ”  -  -  h  <13) 

h  being  the  total  thickness  of  the  shell 
h^  being  the  thickness  of  the  liner 
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For  a  unidirectional  filamentary  layer  which  is  in  a  state  of  uniaxial 
tension,  the  linear  strain  takes  the  following  form  by  the  reduction  of 
the  Eqs.  (A-9)  with  o)  =  (Ref.  6) 
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The  axial  filament  stress 
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The  mean  stress  components  cr,,  a^,  ,  will  be  obtained  by  combining  the 

contributions  of  the  force  per  unit  width  of  the  i-th  filament  layer  (Figure  3), 
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with  those  of  the  matrix/ and  liner,  and  by  summing  over  n  filament  layers. 
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Figure  3.  Forces  in  the  i-th  Filament  Layer  Pair 


where  the  stress  components  in  the  matrix  and  the  liner  are  given  by 
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Substituting  Eq.  (15)  and  (17)  into  (16)  yields 
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Comparing  Eq.  (18)  with  Eq.  (2)  and  associating  6  and  (f)  with  x  and  y 
respectively 
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With  C  determined,  S. .  can  be  calculated  from  Eq.  (3).  The  stress 

ij  13 

in  the  filaments  of  the  i-th  layer  is  obtained  by  substituting  Eq.  (1)  xnto 
Eq.  (15) 
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i  ^ 


Rewriting  in  terms  of  the  stress  resultants. 
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E-  ( 

(7  =-^ 

f.  h  I  d 
1  ' 


2  2 

S.,cos  a.  +S._sin  a.  +  S.  „cos  o.  sina. 
11  i  12  1  13  1  1 


+  N' 


m 

0 


2  2 

S.„cos  O'.  +S„„sin  a.  +  S„  cos  a.  sinQ-. 

1.  ].  X  X 


(21) 


+  N' 


m 

0<l) 


2  2 

S  cos  a.  +S„_sin  a.  +  S  cos  or.  sin  a. 

Xo  X  <u<j  X  OO  X  X 


Similarly,  substituting  Eq.  (l)into  Eq.  (17),  the  components  of  stress 
in  the  matrix  are  obtained 


e<^ 


E 


h(: 

.  2  \ 

[  -  V  ) 

\ 

m 

®23*>'n,^3> 


ii  m  X2  m  22  13  m  23 


a  -  V  ) 

m 


(1  -  1/  ) 
m 


(i  -  ) 

m 


13 


23 


33 


N. 


N  , 


N' 


60 


(22) 


Likewise,  the  components  of  stress  in  the  liner  could  be  obtained  by  replacing 
the  subscript  "m"  with  in  Eq.  (22), 


For  symmetric  filament  pattern:  S  =  S  =  0 

X  o  ^  o 

For  symmetric  filament  pattern  and  symmetric  loading:  ~  ^23 


0 
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Limitations  of  the  Membrane  Theory 


The  membrane  approximation  to  the  particular  solution  of  the  shell 
equations  is  generally  invalid  when  dealing  with  shallow  shells  which  are 
non-isotropic  and  non-spherical  in  the  vicinity  of  the  apex.  Prof.  C.  R. 
Steele  of  Stanford  University  in  his  recent  work  on  orthotropic  shells  of 
revolution  (Ref.  4)  obtained  a  particular  solution,  in  terms  of  Lommel's  and 
Bessel's  functions,  which  is  equivalent  to  the  well-known  "membrane" 
solution  in  the  steep  portion  of  the  shell  but  in  the  shallow  portion  gives 
significant  bending  stresses.  These  bending  stresses  are  not  due  to  edge 
effect  solutions,  but  exist  simply  because  of  orthotropy  and  the  shape  of  the 
meridian  curve. 

The  bending  stresses  are  confined  to  the  immediate  vicinity  of  the 
apex  for  thin  shells;  however,  for  thicker  shells  may  extend  well  down 
along  the  meridian  from  the  apex. 

Expressions  for  the  stress  resultants,  strains  and  deformations  due 
to  the  particular  solutions  are  as  follows  (Figure  4), 


m 


— '‘—z  ver  X 

jU  -  1  n,  V 

o 


for  jj.  I  or  /i  5  with  j  =  1,  2  (23a) 


3  P  r 

N  =  N  +  - — : — j-  ver  x 

01  01  4  sin  0  1,1/ 


for  ju  =  1,  j  =  1  (23b) 


N 


03 


^T.m  P  cot  0  ^ 

- -  ver,  X 

03  47rs  1,  V 


for  M  =  1,  j  =  3 


N 


2a  /  * 

p  rcos  0  /ver^  x 
N,..  + - : — : \ - J - 1/  for  |u 


01  01  sin0 


=5,  j  =  1 


(23c) 


(23d) 


>K 


N 


02 


(ver^  X 

5,  V  , 

- 4 - ■"''*''5,  ' 


for  Id  =  5,  j  =  2  (23e) 


N,.  =  )cos^0  ver  x  +  sin^  0  (  (all  cases) 
0J  01  )  a.  Id.,  V  a  j 

J 


(24) 
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r  cos  6 

- -  n“  vei  X 

,2^/r  “r'' 


(all  cases) 


in  which  6  is  the  vertex  angle  for  b  =  1  (for  b  1,  -  0). 

a 


The  stress  couples  are 


hcos<^i  f  b  *  / 

M„.  = - ^N^(C_cos(^.  vei  x+C,„fvei'  x[  (all  cases) 


hcos(/i  f  b  *  ) 

M,.  = - ^  nT!  iC,^cos  (/>  vei  x+C„„— vei  x[  (all  cases) 

03  03(  12  a  p.,  1/  22  2  n.,  u  j 


where 


03  •  2  , 


3  =  1,  2,  3 


= 

^3 


P^2 - 2  j  "  ^ 

j  r^sin  0 


r^sin  <!> 


3  =  2,  3 


=  A  Ic  -  C  n”" 

ej  hA  \^22^0j  ^12  0] 


3  '  1.  2,  3 


p  r  cos  <j>ds  ,  T_ 
r  z 


p  rsin0(is,  = - 

S  3  27r 
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The  argument  of  the  Lommel  functions  is  given  by 


where 


A 


S1S2 


-  C 


2 

12 


and 


000 


(31) 


(32) 


(33) 


are  the  normalized  distances  with  R  =  max  r 

o 


IX. 


6  -  b 


;  ^ 


4  +  b 


2  b  ’ 


- 


2  -  b 


2  4  11 


V 


Z  =  D  r^^  |l  +  0(1^ ^)j  for  b  =  1  meridian  curve 


(34) 


(35) 


(36) 


The  value 

b  =  1  gives  the  class  of  shells  that  are  conical  at  the  apex; 
b  =  2  gives  the  class  of  shells  that  are  spherical  at  the  apex; 

and 

b  >  2  gives  the  class  of  shells  with  zero  curvature  at  the  apex. 
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The  Lommel  functions  (Ref.  7)  which  are  analogous  to  Thomson 
functions,  expressed  in  terms  of  real  and  imaginary  parts  are  as  follows: 


(37) 


From  the  asymptotic  behavior  of  V 
for  fixed  S  while  k-*-oo,  the  result 


(A)  and  V*  (A)  for  large  A  one  has, 

V  fJ.,V 


^  £ 


m 


X.  ~M  ,  .~M  ~0 

J 


(38) 


The  behavior  of  the  Lommel  functions  for  /i  =  1  and  2  is  shown  in 
Figures  5  and  6. 

The  results  are  uniformly  valid  throughout  the  shell  and  reduce  to 
the  shallow  shell  solutions  for  small  A  and  approach  non-shallow  solutions 
asymptotically  for  A  »  1. 


Application 

The  results  of  the  above  formulation  as  applied  to  a  shell  of  revolution 
with  two  clamped  edges  at  S  and  S  -  S^(S^  rigid  rings  can  be 

advantageously  represented  in  nondimensional  curves  for  maximum  stresses 

at  the  edges  depending  only  on  v,  b,  C22/C22  parameters.  These  curves 
represent  the  behavior  of  shells  of  wide  geometric  variation.  The  curves 
presented  for  the  various  numerical  values  of  the  above  parameters  will 
enable  the  stress  analyst  to  interpolate  for  his  given  shell  configuration. 
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8  0 


For  the  case  of  a  spherical  shell  subjected  to  uniform  pressure,  p^. 


b  =  2,  M  =  2, 


The  Figure  7  taken  from  Steele's  paper  (Ref.  4) 


shows  the  nondimensional  curves  for  ^22  ”  1/3,  and  7/8  with  u 

representing  the  family  of  curves.  Curves  with  the  M  represent  bending 
stresses,  while  those  with  N  represent  direct  stresses. 


The  behavior  of  the  filament-wound  shell  can  be  approximated  by  taking 
local  values  of  x  and  u  along  the  meridian  of  the  shell,  Hartung  (Ref.  8). 

This  is  equivalent  to  treating  the  shell  as  if  it  were  composed  of  many  seg¬ 
ments  each  of  which  has  constant  x  and  u,  and  then  requiring  the  stress  be 
continuous  from  one  segment  to  the  next.  The  x  and  u,  which  are  necessary 
to  compute  the  Lommel  functions  are  determined  from  expressions  (31)  &  (35). 
The  value  of  the  argument  of  the  Lommel  functions  depends  upon  the  fila¬ 
ment  size  and  spacing,  the  moduli  of  elasticity  of  the  filament,  matrix  and 
liner  as  well  as  the  applied  loading. 

With  the  argument  of  the  Lommel  functions  x,  determined  from 
expression  (31);  the  modified  and  N,  can  be  found  from  (23)  and  (24). 

These  are  then  used  with  (21)  and  (22)  to  calculate  the  modified  "membrane" 
solutions.  The  bending  stresses  produced  because  of  orthotropy  and  the 
shape  of  the  meridian  curve  can  be  calculated  by  substituting  the  values  of 
the  bending  moments  obtained  from  the  equation  (26)  into  equations  (52) 
and  (53). 


Determination  of  the  Winding  Angle 


The  filaments  are  deposited  by  planar  winding  (Figure  8). 

The  winding  angle  /3  between  the  meridian  and  the  filament  is  found 
from  the  definition 

cosB  =  t  'T  -  t  T  +  t  T  +  t  T  (39) 

XX  y  y  z  z 

where  t  and  T  are  unit  vectors  tangent  to  the  meridian  and  the  filament 
curves,  respectively.  The  components  of  these  vectors  can  be  obtained 
from  the  following  expressions  with  primes  denoting  differentiation  with 
respect  to  z. 
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Fig,  7a  Stress  Resultants  at  Edge  of  Pressurized  Dome  Clamped  to  Rigid  Ring 


Stress  Resultants  at  Edge  of  Pressurized  Dome  Clamped 


Figure  8.  Geometry  of  Planar-Wound 
Shell  of  Revolution 


r'  cos  6 


t  =  cos  { t,  x) 

X 


-  cos(t,  y) 


t  =  COS  (t,  z) 
z 


X' 


Vx'^  +  y'^  +  1 


x'^  +  y'^  +1 


i. 


T  =  cos(T,  x) 


T  =  cos(T,  y)  = 


’  T  =  cos  (  T,  z)  = 


x'^  +  y'^  +  1 


i 


1  +  r' 
r'  sin  0 


(40) 


1  +  r'‘ 
1 


1  +  r' 


The  projection  of  the  filament  curve  on  the  (y,  z)  plane  is  given  by 


2  2 
X  =  -^ r  -  (y^  +  z  tanT) 


(41) 


where  y  is  the  offset  of  the  winding  plane.  R.  Hartung  (Ref.  9)  obtained 
a  gener^  expression  for  the  winding  angle  jS  as  follows 


tan 


2 


/3  = 


|rr'  -  (y^  +  z  tan  7)  tan  7 

^22 
+  (tan  7  “  ) 

|r^  -  (y^  +  z  tan  y)^! 

(1  +  r'^) 

2 

r 

"2 

-  (y  +  z  tan  7) 

^  e  '  J 

(42) 


For  the  special  case  of  a  spherical  shell  (Figure  9)  wound  along  its 
geodesics,  by  substituting  the  following  expressions 

r  =  -  z^ 

r  =  Rcosi// 
z  =  Rsin^// 

also  y  *  0 
•^e 


r'  =  ^  =  -  —  =  -  tan^  (43) 

dz  r 


.34 


sin  /3 


sin  y 
s  in  (j) 


Figure  9.  Planar-Filament-Wound  Spherical  Shell 
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the  equation  (42)  becomes 


2,  rRcos(//(- tarn//)  -  (R  sin  (i'tanY)  tan 
tan  p  -  - -  2  r  2  2 

i/yH  R  nofi  f// 


9  r  9 


which  reduces  to  a  simple  relationship  among  the  three  spatial  angles  /3, 
<l>>  7 


sinjS  = 


sin  y 
sin  0 


where  0  =  ^ 


(43) 


Bending  Stresses 


Assuming  that  the  shell  thickness  is  made  up  from  a  large  number  of 
layers  of  the  same  elastic  properties  and  further  assumirg  that  the  variation 
of  the  moduli  of  compliance  S. .  along  the  meridian  of  the  shell  is  slow,  for 

unsymmetric  loading  (Ref,  6), 


e  =  e  =  Ik. 

yd  <P 

^xy  ^6<l>  ^  9^ 


(44) 


where  the  changes  in  the  meridional  and  the  circumferential  curvatures  and 
the  twist 


ax^ 

'"(t>  "  as 


K 

9 


I  /  9Xq 
r  \  a0 


60  r j  9  0 \  r  /  r  96 


(45) 
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the  stresses  are  obtained  from  Equation  (2) 


"9 


(46) 


The  stress  couples  (Figure  10) 

•h/2 
-h/2 


M^,  M^, 


rh/2 


(47) 


yield 


“a 

,3 

•^ii 

<=12 

<=13 

% 

h 

"  I2 

=12 

S2 

S3 

% 

(48) 

“ao_ 

.'"is 

^23 

S3_ 

In  symbolic  form 


[M]  =  [C]  [K] 


(48a) 


M  =  l^\  [C]'^  [M] 


(49) 


•1  . 


where  [C]  is  the  inverse  of  lC]  and  is  given  by 
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Axis  of 
Revolution 


Figure  10.  Stress  Couples  Acting  on 
Differential  Element 
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[C] 


-1 


r22^33  ^23 


)  (< 


-  C,„C  ) 

(c.  c 

-  c  c 

13  23 

12  33/ 

\  12  23 

13  22 

-c\\ 

(c,„c. 

-  c. 

11  33 

13/ 

\  12  13 

11  23 

(c, 

-c?  ^ 

V  11  22 

12/ 

where 


A  = 
c 


^12  ^13 


^2  S2  S3 


<=13  S3  C33 


(50a) 


Using  expression  (3)  in  Equation  (49)  and  substituting  in  turn  in 
Equation  (44) 


e 

0 

■  ‘&) 

\h  / 

1 

_ 1 

S. 

S. 

S. 


s, 

S. 


s, 

s, 

I 

s. 


1 - 

CD 

1 _ 

(51) 


Substituting  in  Equation  (15),  the  axial  filament  stress  in  obtained. 


“  M  “5"  1  1 +S,„sin  a.  +S._coso'.  sincv. 

I.  i  12  1  13  1 

1  \h 


+  M,(s.  cos^o.  +S  sin^o.  +  S .  „  cos  c.  sincv.  )  (52) 

0\  12  1  22  1  13  11/ 


(2  2 
S  cos  a.  +S„-sin  a.  +  S„„  cos  a.  sina. 

Xo  1  2o  1  1  1 


(50) 
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Similarly^  Equation  (17)  yields  the  stresses  for  the  matrix 


O', 


m  \h  /\1  -  V 


E  \  , 
m  \  I  _ 


M.  S._  +  1/  S 
6\  12  m 


u)"^( 


S  +1/  s 

22  m  12 


m 


+  +,/  S 


14 


23  m  13 
E 


m 


m  \h V  \  1  -  i. 


2  I  [^0 
m , 


(S  +r  +1/ 

\  11  m  12/  (j)\  12  m  22y 


(53) 


O'/ 


04* 


1/12 


m  \h 


3/2(1 


E  / 
-  gi-  ■  / 

+  v  )\ 

m 


For  symmetric  ■winding  and  unsymmetric  loading:  S  =  S  „  =  0 

1 0  ^  o 

For  symmetric  winding  and  symmetric  loading:  =  S  _  =  M  ,  =  0 

13  23  Q(p 


Similarly,  the  components  of  stress  in  the  liner  are  obtained  by  replacing 
the  subscript  "m"  with  in  Equation  (53). 


The  bending  stresses  can  thus  be  computed  in  terms  of  the  stress 
couples  Mg,  and  the  twisting  moment  applied  on  the  edges  of  a 

shell  element.  For  the  edge  effect  behavior  of  a  shell,  the  twisting 
moment  plays  a  secondary  role  and  may  be  neglected;  however,  the 

stress  couples  play  an  essential  role.  The  principal  one  of  them  is  M 

0 

The  stress  couples  will  be  furnished  in  the  section  "Edge  Effect  Solutions  " 


in  terms  of  the  membrane  stress  resultants,  N™  and 

0  0 
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Reduction  of  the  Equilibrium  Equations 
to  Two  Differential  Equations  of  the  Second  Order 

Denoting  the  angle  of  rotation  of  a  tangent  to  a  meridian  by  X 


dw  , 
X  =  + 

ds 


u  1  /  dw  ,  \ 

r^  =  r^  (■  d</.  i 


(54) 


and  introducing  U  =  ^2^(5  ’  considerable  simplification  can  be  achieved  in 
the  following  formulationTRef.lO). 

In  the  absence  of  surface  loading,  the  equilibrium  equations  become, 

N  ,  =  Q  ,  cot  (/) 

0  0 


,  d(rQ  ) 

r  N/,  +  r  N  =  — = - ^ 

1  6>  ^2  (f,  sin0  d0 


(55) 


N 


N. 


r. 


-U  cot  0 

dU 

d0 


(56) 


For  axisymmetrical  loads,  with  the  substitution  of  the  strain  components 
of  the  middle  surface  in  terms  of  the  deformations. 


J_  /  ^ 
r^  r  d0)  ^ 


(w  +  u  cot  0) 


(57) 


the  stress  resultants  and  stress  couples  become. 


N.  =  h 


N  ,  =  h 


'll 


(w  +  U  cot  0)  + 


^12  /  ^  du 

-  w  +  XT 

Tj  \  d0 


)] 


rr*w  +  licot*)  +  ^(w  +  ^)] 


(58) 
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(59) 


-n 


^(t>  12 


1  ,,  ,  ,  ,  '"12  dx 

C , .  —  X  cot  0  + - 7-7- 

11  d(f 

_  1  ^  ^  S2  dX 

C  — X  cot  0  +  - — 

12  d0 


To  establish  the  first  equation  connecting  X  and  U,  the  Equation  (58) 


gives, 


w  + 


1 


du  ^1 


d0  h  /  2  ^  ^  ' 

\  12  11  22; 


W  +  U  cot  0  = 


K1S2  -  C12) 


2  \  S2  -  ^12, 


(60a) 

(60b) 


Eliminating  w  from  these  equations,  and  using  expressions  (3)  for  the  moduli 
of  compliance  S, 


du  1 

=  h 


^^1^22  "  ^2®12^  ^0  “  ^^2^11  ^1^12^  ^0 


(61) 


Differentiating  Equation  (60b) 


4.  ^ 


u  ^  dw  d 


.  2  ^  d0  d0 
sin  0 


)] 


(62) 


The  derivative  of  u  can  be  readily  eliminated  from  Equations  (61)  and  (62) 
to  yield 


dw  _ 
d0  ^  d0 


h  \®11^0  ®12^0 


)l 


'  ^^2®11  ■  ^1^12^  ^0 


Substituting  expressions  (56)  for  N  ,  and  N  ,  an  equation  relating  U  and  X 

0  0 

is  finally  obtained. 
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®lll 

j  d^U  dU 

fll  d  , 

h  ‘ 

[r,l 

d* 

h  d^  ' 

[rj  d(i) 

d  /I 


cot  (/)  (  a  _  r  S  ) 
hr^  ‘'^2^11  "^1^12' 


+  U 


cot0  ‘'®12  ,  „  H  <i  lA 

- +  d*y 


h  d<f) 


d0 


_  co^,  -  r  S  ) 

hr  ^1^22  2^12' 


•  2  ,  h 

Sin  (j) 


=  -r^X 


(63) 


For  the  isotropic  case: 


E 


1  -  r 


2  ’  ^12  ^ 


V 


1  -  V 


^11  =  ^22  =  E 


;  s 


12 


V 

E 


(64) 


and  hence  the  Equation  (63)  reduces  to  Equation  (315)  of  Timoshenko , (Ref.  10). 

The  second  equation  for  U  and  X  is  obtained  by  substituting  expressions(59) 
for  Mg  and  in  the  following  moment  equilibrium  equation. 


jj(rM^)  -  Mg  cos^  -  rQ^  =  0 


(65) 


which  in  turn  yields 


r  2 

_2  d^  J_ 

2  ^^2  r. 
r dp  1 


+  r. 


dX 

d4> 


cot  ^ 


^/n  ^3^1  £i2  X  [/, 

d'f'VV  '■l ’’d*!  +  0^2  rj  [I  ■ 


2 

cot  (/) 


eo,  i  d  n  )  .  1  ^  3  ^  u  /m 


(66) 


22  'h 


The  Equation  (66)  reduces  to  Equation  (316)  of  Timoshenko  (Ref.  8)  for 
the  case  of  isotropy. 


Thus,  the  problem  of  bending  of  a  shell  of  revolution  under  stress 
resultants  and  couples  along  its  edge,  a  parallel  circle,  is  reduced  to  the 
integration  of  Equations  (63)  and  (66).  Since  the  stresses  produced  by  the 
stress  resultants  and  couples  applied  at  the  edge  of  the  shell,  present 
themselves  as  a  rapidly  decaying  state  of  stress  when  the  angle  (fi  is  not 
too  small,  a  simplified  solution  can  be  obtained  by  means  of  Geckeler 
approximations  applied  to  spherical  shells.  Neglecting  the  functions 
U  an(d  y  anci  their  first  derivatives  in  comparison  with  the  second 
derivatives,  the  Equations  (63)  and  (66)  become. 
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X 


(63a) 


A 


d^X 


2 

h  ^1 


®11  ^2 


12 


d(^)"  ^^2  "^11 

These  equations  can  be  combined  into  a  single  differential  equation. 


1  (^11^22  ~  ^12) 


U 


(66a) 


+  4k^X  =  0 

d04 


=  -2[- 
h  V2> 


Ki^22  -  ^12) 


11 


(67) 


(68) 


Here  k  is  a  function  of  0  angle.  A  satisfactory  approximate  solution  of 
Equation  (67)  can  be  obtained  by  replacing  k(0)  by  a  certain  constant  average 
value  as  suggested  by  Timoshenko  (Ref.  10). 
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Edge  Effect  Solutions 


Every  state  of  stress  arising  in  a  shell  as  a  consequence  of  the  fact 
that  it  has  a  line  of  distortion  is  called  an  edge  effect.  The  lines  of 
distortion  may  be:  1)  the  edges  of  the  shell;  2)  lines  along  which  dis¬ 
continuities  of  the  components  of  the  external  surface  loads  occur; 

3)  lines  along  which  the  curvature  of  the  middle  surface  changes  abruptly; 

4)  lines  along  which  the  rigidity  of  the  shell  or  its  thickness  undergoes 
sudden  changes. 

A  complete  study  of  the  state  of  stress  consists  of  the  membrane 
analysis  of  the  shell  and  a  determination  of  the  local  stresses  occurring 
near  lines  of  distortions  due  to  edge  effects.  The  membrane  analysis  gives 
an  approximately  correct  picture  of  the  state  of  stress  and  deformation  of 
the  shell  only  at  sufficient  distances  away  from  the  lines  of  distortion. 

The  edge  effect  presents  itself  as  a  rapidly  decaying  state  of  stress  away 
from  the  edge  of  the  shell. 

A.  Gol'denveizer 's  Solution 

Assuming  a  system  of  curvilinear  coordinates  7^,  and  7^  has  been 
constructed  on  the  shell  (Figure  11),  in  which  the  line  of  distortion 


7^q  =  constant 


Gol'denveizer  obtained  (Ref.  1). 


1  d^X  ,  12(1  -  v^) 


4  4 

A  d7j 


0 


^2 


(69) 


The  quantities  r'  and  v*  represent  the  radii  of  curvature  of  the  normal 

1  Cl 

sections,  cut  along  the  coordinate  lines. 


First  fundamental  form  of  the  theory  of  surfaces  is 
(ds)^  =  A^(d7^)^  +  2ABCOS  d7  jd72  +  B^(d72)^ 


(70) 
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Figure  11.  Differential  Surface  Element 
in  Curvilinear  Coordinates 


where 


0,7 

~~A~ 


0,7, 


=  cos  Q ; 


Sr" 

0 

;  B  = 

0 

57  1 

^72 

(71) 


Tq  is  the  position  vector,  ^q('Vj^,  7 2^ 

Q  is  the  angle  between  the  coordinate  curves. 

For  orthogonal  coordinates, associating  7  with  (/;  and  7 2  with  6  ,  the 

surface  is  referred  to  lines  of  curvature  coordinates,  since  the  lines  of 
curvature  of  a  surface  of  revolution  are  its  naeridians  and  parallels.  Hence, 

the  quantities  r'j  and  r'  become  r^  and  r^,  i.e.,  the  principal  radii  of 

curvature,  and 

A  =  r^,  B  =  r,  0=0  (71a) 

Then, 

(ds)^  =  r^(d(^i)^  +  r^(d0)^  (70a) 

Comparing  the  Equation  (67)  with  Equation  (69) 


S1A2 


-  S 


12 


'11 


and 


(72) 
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Therefore,  the  stresses  and  deformations  in  orthotropic  shells  can  be 
calculated  from  the  Gol*denveizer‘s  edge  effect  solutions,  within  its  range 
of  validity,  by  means  of  the  above  relationships. 

The  membrane  state  of  stress  is  determined  by  the  quantities: 


N' 


m 
21  ' 


m 


m 

V 


m 

w 


X 


m 

1 


and  the  edge  effect  by  the  quantities: 

sic  5}c  s|c  s|c 

Mj,  W  ,  Xj,  Qj 

Hence,  the  stress  resultants  and  displacements  of  the  complete  state  of 
stress  may  be  written  in  the  form: 


iNf 

vl>- 

'I'* 

=  ; 

u  = 

m 

u 

>“2 

❖ 

=  ; 

V  = 

m 

V 

“l2 

“l2 

=  M21  =  0  ; 

w  = 

m  .  * 

w  +  w 

(73) 

“21 

*^1 

vt. 

=  q; 

h  = 

ym  ,  y* 

Xi  +  x^ 

^2  =  ° 


The  formulas  for  the  edge  effect  theory  are  outlined  below.  Here,  and  in 
the  following  pages,  the  subscript  0  indicates  the  edge  (Figure  12). 


48 


1)  For  the  lower  edge,  y  ^  s  7 


Ehw'  =  cos  kg  (7  j  -  7  jq)  +  /"2  sinkg  (7  ^  -  7  jq)]  1  "'^10^ 


3(1  -  v^) 


7^{[^  1  +r2)coskg(7j  -7j(,)  -  [Tj  +r2]sinkg(7j  -  7  j„)} 


>|c  1  * 

M,  =  -M^ 
1  V  2 


-^j|^=S=^[r2COSkg(7^  -7j„>  -r^sinkg(7j  -7io) 

2 


*7  1^  1  -■r2jcoskg(7j  -7jo)  -  |Oj  +  r 


'1  4 


12(1  -  ^2 


X  sinkg(7  -  7  jq)|  1  "^10^ 


^2  =  -  ^  |r^  coskg(7  j  -  7^q)  +  ^2  ■  '^10^  1  ”^10^ 


2)  For  the  upper  edge,  y^  ^ 


Ehw"  =  jr^  cos  kg  (7^  ■  ”^10^  '*’  ^4  sin  kg  (7^  -  7^^)]  e  l  ”^10^ 


[^3  ■  ■r4j«>skg(7i  -  7jo)  -  [rj  +  rj 

y^''2 


X  sin  kg  (7  ■  7  ^q)|  e 


-kg(7i-7io) 
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>}=  1  * 

M,  =  - 

1  v  2 


2r2V3(l  -  v‘^) 


[-  cos  kg  (7^  "  r^sinkg(7^  -  T^q)] 


X  e 


-kg(Ti-Tj(,) 


?ic 

«1  '  4 


* 


N. 


Vl2(l-r2) 

X  sinkg(7j  ■'>'10^}®  ^10^ 

=  -  —  [rgCOskg(T^  ■  '''10^  '*’  ^4  sin  kg  (7  j  -  'Vjq)]  e  l  10^ 

«/  I 


(75) 


where  r^.  F^  are  arbitrary  functions  of  72 


and 


kg  = 


3(1  -  r“) 


2.}'^  A 


20 


(76) 


A  and  r'  are  the  values  of  the  quantities  A  and  r'  at  the  edge  of  the  shell. 

0  <y  U  ^ 

The  Edge  Effect  Near  a  Clamped  Edge 

1)  Upper  Edge  (7^  ^  7^^  or  (j>  ^  <l>^) 


The  boundary  conditions  at  the  supported  edge,  =  7  are 
u=0;  v=0;  w=0;  X^=0  (77) 


The  first  two  of  these  four  conditions  are  tangential,  the  last  two  are  non- 
tangential.  The  tangential  boundary  conditions  have  been  fulfilled  by  the 
membrane  theory,  since  by  Expressions  (73) 


ni  n 

“0  "  “0  '  “ 


^0  =  ° 


(78) 
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The  non-tangential  boundary  conditions  may  be  written  in  the  form. 


m 

Ehw^  =  EhWp  +  Ehw^  =  0 


EhX^O  ^  ^ 


For  7  j  =  7  using  Equation  (75) 


EhWo  =  Tj 


EhX 


V 


10 


3(1  -  v") 


r  -  r 

3  4 


Hence 


^3  = 


-  Ehw 


m 

0 


^4  = 


Ehw“  + 


0  4^ 


EhX 


m 

10 


3(1  - 


(79) 


(79a) 


Obviously,  for  the  membrane  analysis  of  the  shell,  w  and  1^20 ^1 

quantities  of  the  same  order,  and  on  the  right-hand  side  of  the  above 
expression,  therefore,  the  principal  part  will  be  played  by  the  first  term, 
so  that  we  may  assume  as  an  approximation  that 

r  =  -  Ehw“  (7  9b) 

4  0 

Then,  for  the  computation  of  the  stress  resultants  and  moments  of  the  edge 
effect  due  to  loads  of  rotational  symmetry 
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=  -  -y  [r^  cos  kg  (Tj  -  7jq)  +  r^sinkg  (7^  -  7^Q)]e  1  10^ 


""  2r'V3(l  - 


[-  cos  kg  (7  ■  '^10^  ^3  si"  kg  (7  ^  -  7^q)] 


A  and  r'  are  the  values  of  the  quantities  A  and  r'  on  the  boundary  of 
0  U  ^ 

the  edge. 


Substituting  for  and  F^  in  Equation  (80) 

=  -  —  |-  Ehw^coskg(7  j  -  7  jq)  -  Ehw^sinkg(7j  -  7j^Q)je  1  ”^10^ 


cos  kg  (7  j  -  7  ^q)  +  sin  kg  (7  10^  ®  ^ 


°  ^  '  '  2r'V3(‘l  -  vh  '  ^10>  '  Ehw^^l"kg(T^  -  Tjp) 


^^-kg(Ti-Yi„) 


h  ^‘^'"0 
2t'  Vsd  -  v^) 

<u 


|coskg(7j^  "  ~  si"  kg  (7  -  7j^Q)]e  ”^10^ 


To  eliminate  observe  (Ref.  1) 


1  ^  ^  1  >  -D 

^  V  1  d  B  m  w _ 

B  d  72  -^B  d  7  2 

2 


1  /,,m  Tyrm) 

EhK  -■'“l 
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Since  at  t  ^  =  T 

m  m  ^ 
u  =  V  =0 


w. 


m 


20 

Eh 


(“20  -  "“loj 


Thus,  Equation  (81)  may  be  written 


N. 


-  -^(n“  -  |coskg(Tj  -  7jq)  +  sinkg(T^  "  '^lO^j 


r'  (n^  - 

h  20\  20  10/ 


>{= 

^  2 

M,  =  -  ^  I - ^ 

^  ""  2r'  Vsd  -  v^) 


|coskg(T^  ■  '^lO^  '  -  10>1 


(82) 


^e-kg(Ti-7 


10 


) 


In  the  above  formulas,  the  quantity 


should  be  set  equal  to  unity. 


because  such  an  assumption  had  been  adopted  by  Gol'denveizer  in  the 
general  formulation  of  the  theory. 


By  choosing  7^=0  and  72  =  ®  »  'the  surface  is  referred  to  lines 

of  curvature  coordinates,  since  the  lines  of  curvature  of  a  surface  of 
revolution  are  its  meridians  and  parallels. 


Hence,  the  quantities  r'^  and  r^  become  r^  and  r^,  i.e. ,  the 
principal  radii  of  curvature. 


Rewriting  Equation  (82)  in  the  notation  of  the  membrane  theory  presented 
earlier  in  this  report,  the  following  expression  Equation  (83)  for  the  upper 
edge,  0  =*^0’  obtained. 
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N, 


-T.m 


|cosX((^  -  +  sinX((/)  -  '/'q)]© 


(83) 


*  ^0 


“e  - 


<^  =  0 


= 


Wsd  -  v‘ 


—  cos  X{<t)  -</)^)-sinX((/)  -</»q)  e  ^  0^ 


where 


V 


3(1  -  v^)  ^10 


{I 


hr, 


20 


(84) 


2)  Lower  Edge,  (7  S  7^^^  or  S  (f)^) 


* 


N, 


=  -  (n“  -  rN^]  [cosX  (0  -  <l>^)  -  sinX  (0  -  (j) 


e  0 


h 


M. 


Ne  - 


0=0 


M  ,  = 


(fi  V 


2'^3{1  -  v‘ 


—  |cos  X  (0  -  </) q)  +  sin X  (0  -  0q)|  e 


(85) 


X(0-0q) 


For  spherical  shells,  the  position  vector, 

r„  =  R  sin  6  cos  d  e  +  R  sin  0  sin  0  e  +  R  cos  0  e 
0  ~x  ~y  ~z 

where  e  ,  e  ,  e  are  unit  vectors  along  the  coordinate  axes  x,  y, 
~x'  ~y'  ~z  ^ 


z,  and 
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Associating  7^  with  <j)  and  'S ^  with  Q 


,  R  cos  0  cos  6  e  +  R  cos  0  sin  S  e  -  R  sin  0  e 
0rjS0  ^  ~y  ~z 


A  = 


S0 


=  R 


and  similarly  B  =  r  =  Rsin0 


Since  r  =  r  *  R 
1  ^ 


(84a) 


The  Edge  Effect  Near  a  Hinged.  Supported  Edge 
1)  Upper  Edge  (7^  ^  7^^  or  0  ^  0^) 

The  boundary  conditions: 

=  0  ;  V  =0  (tangential) 

w  =  0  ;  =  0  (non-tangential)  (86) 

m 

Ehw^  =  Ehw^  +  EhWp  =  0 
* 

=  0 

From  Equation  (75),  at  7  =  7^^ 

r^(l^  =  -Ehw“  ;  r^(y^  =  0 

Substituting  in  Equation  (80)  and  changing  the  notation. 
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0  r„  0 


M 


* 


1  * 
-M^ 

V  e 


2  m 
Eh  Wq 


2r^^30~^‘ 


-X (0-0)  . 

e  ^  ^0  sinX(0 


where  X  is  given  by  expression  (84a) 


(87) 


2)  Lower  Edge,  or  </•  5  0^) 


THU 

.  Ehw^  ^  . 

=  - ^y^e^^'^"'^0^cos  X(0  -  0q) 


X:  1  * 

M  ,  =  -  M„ 

0  V  0 


-rp,  2  m 

*0  Mt-iJ . 

- 1  .  ■  e  0  sinX(0  -  0„) 

2rW3(l  -  v^) 


where  X  is  given  by  Equation  (84a). 


(88) 


B.  Asymptotic  Solution 

The  equilibrium  equations  for  the  case  of  rotational  symmetry  (Ref.  2), 

2 


Eh"^ 

—  1 

fr^l 

l-l 

L.  sin0 

12(1  -  v^) 

ds  ' 

1  dsj 

1/ 

^  ^1 

rH  sin  0  =  -  rV  cos  0 


1 

Eh 


1(4  H  -  f 


(89) 


+  xsini,  =  ^ 


+  cos  0  (Vsin0  -  I'rp^) 

id 


where,  X  is  the  rotation  of  the  normal  to  the  shell  and  V  and  H  are  the 
axial  and  radial  components  of  the  stress  resultants  (Figure  13), 
respectively. 
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The  axial  and  radial  components  of  the  stress  resultants. 


V  =  N  ,  sin  <j)  -  Q,  cos  (j) 

0  0 

H  =  N  ,  cos  0  +  Q  ,  sin  0 

0  0 


(90) 


Similarly,  the  deflection  of  a  point  on  the  midsurface  can  be  given  by  the 
radial  component  6  and  the  axial  component  v. 


6  =  w  sin  0  +  u  cos  0 

V  =  -WCOS0  +  usin0 

For  homogeneous  equations. 


(91) 


P  =Pj.=P  =Ptt=0 
'^<j)  H 


Introducing  a  non-dimensional  quantity. 


Eh 


the  homogeneous  equations  are  transformed  into  the  form  (Ref.  2) 


d  /  dX\  /  sin  0  ^  cos  X  -  ~ 

ds  \  ds  /  \  r  J  h 


^sin  0  =  0 


^  H-  (z.  xsin  0  =  0 


(89a) 


Assuming  a  solution  in  the  form  of  a  series  expansion 

;^  =  ^  +  h^  +  — -j 


X  =  [  Xq  +  h^'^X^  +  hXg  +  — -J 


_  €h 


(92) 


where  ?  ,  X  .,7^.  are  bounded  functions  of  s  that  are  completely  independent 
3  3 

of  h. 
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Substituting  into  the  homogeneous  equations  and  equating  the  co¬ 
efficients  of  h's  equal  to  zero  and  retaining  the  leading  terms  in  the 
expansion,  the  solution  becomes 


X  =  ReCe*'^  +o((h/r)‘^^) 

where  C  is  a  complex  constant 


(93) 


3(1  -  i'^) 


1/4 


-  ^0>] 
V  ‘’”’20  I 


(94) 


'■20  °  ‘'■2>edge 


T  r 

1  2 


For  thin  shell,  r^,  r^  »  h,  say  —  ,  ^20  for  the  simple  edge 

effect  solution  to  be  valid 


12(1  -  v^)r^sm% 


1/4 


cos  </) 


^5 


(95) 


For  sphere,  r  =  Rsin</»,  —  =  1 

^1 


12(1  -  v‘^)Y? 


h‘ 


1/4 

tan  (/)  ^10 


(96a) 


For  pressure  vessels,  simple  exponential  "steep  shell"  solution  gives 
valid  results  for 


12(1 


2vo2 

i'  )R 


1/4 


tan  0^3 


(96b) 
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Stress  Resultants  and  Displacements  for  Upper  Edge 


X  =  (CjCOsI  +  CgSinDe 


(93a) 


M*  =  - 

*  12(1  -  v^) 


+  „xS2S± 

ds  r 


(97) 


Eh  ,,  dX 

12(1  -  '  05 


(97a) 


where 


»/  ^  ^ 
ds 


>}:  >!« 

M.  =  VM, 
t)  0 


In  the  absence  of  surface  loads,  the  horizontal  equilibrium  is. 


N 


6  ds 


(rH) 


(98) 


With 


rH  = 


Eh 


N 


* 


_  d^ 


Eh  d  X 


®  4r^^'  d?  ^ 


(98a) 


Since 


N*  =  N*0(N/h7?)  ; 


Eh 


(where  O  denotes  order  of  magnitude) 
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the  Equations  (97a)  and  (98a)  become 


12(1  - 


[(C^  -  C^)cos|  -  (Cj  +  C2)sin§]e‘ 


^  Eh  r  -S 

Ng  =  _££_  |(c^  +  C2)cos5  +  (Cj  -  Cj)sin5|  e 

^''2^ 


Membrane  Deformations 


Horizontal  deformation: 


'  -L-LJ.  ^  /  , 

)  1  L-xn 

~  ^  d  ""  Eh  re  ■ 


(100) 


for  uniformly  distributed  loading,  p 


m  pr  r  r 

=  __2  2  ._2  .  i; 

r  2Eh  r. 


pr  cos  </>  *  r  d  /  2 


oTT^u  •  2  ,  )cos  <t>  ds  \  r. 
2Ehsin  (f)  I  '  1 


(101) 


Boundary  Conditions  for  Fixed  Support 


m  * 
6+6  =0 


m  ^ 

X  +  X  =  0 


yields 


(102) 


1 


2r2?' 


e"^[(C^  +  C2)cos€  +  (C2  -  C^)sin?] 


X  =  (Cjcos?  +  C2sin?)e 


(102a) 
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At  the  support 


*  ^  m 

—  =  - - - (C  +  C  )  =  -  — 

■^0  ‘  ^  -0 


^  m 

>‘0  =  s  =  -*0 


Hence, 


o  i 

0  ^20  m 

r  ^  0 

0 


Substituting  in  the  Equation  (99) 
„  3  r/ 


12(1  -  v^) 


72  6  ^r-- 

Eh 

(  0  20 

2r^? 

A  0 

<^20.,  .  ..m\  .  ,  /^C^20.A  .  J  ■ 

^0  V  V  -^0  /  J 


^*0  “'20  ,/  ^  ■  el  -I 

- C  +  2X_  Jsin§  e 

*•0  “/  j 


(103) 


?  I  cos?  + 


Since  for  membrane  theory  6^  and  ^20^^  same  order  of 

>{c  ^ 

magnitude,  the  expressions  for  and  can  be  reduced  to  the  following 
for  general  loading  of  a  shell  of  revolution, 

3  2/6?^r„A 


Eh^  _/.^/^0^2o\p  ^  .  ^-1  -€ 

-  )  I - )[cos?  -sin?Je 

6(1  -  \  ^0  ' 


6 

Eh  0  ^20  r  _  ,  .  • 

- cos?  +  sin?Je 

2  0 


(103a) 
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Rewriting, 


M 


hN^-rN^ 

Mg  ^  ®  ? 


(j)  V 


2V3(1  ~  i'^) 


-  |cos(?  -  ^  )  -  sin(?  -  ^Q)fe 


(104) 


N. 


* 


[  0  0. 


-(?-r  J 


|cos(?  -  ?q)  +  sin(?  -  ?Q)}e  0 


20 

Gol'denveizer  had  adopted  the  assumption:  -  ~  1  in  his  formulation. 

^2 

Hence  a  comparison  of  the  Gol'denveizer  and  asymtotic  solutions  reveals  that 
the  two  methods  for  "Edge  Effect”  theory  yield  identical  results. 

For  spherical  shells. 


9  1/4 

?  =  ?  -  =  [3(1  -  )]  (0  -  0q)^j^ 


(94a) 


Deformations  for  a  Shell  of  Revolution 
Under  Symmetric  Loading 


Membrane  Deformations 


1  /du  .  \ 

(d0  7 

=  ^U  cot  0  +  w^ 


(105) 


(57) 


Equating  Equations  (57)  and  (105) 


—  (ucot0  +  w)  =  S^^Og  + 
2 


(106a) 


'  ®12“9 


(107b) 
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Eliminating  w,  the  differential  equation  for  u  in  terms  of  stress  resultants  is 

-  “cot#  =hK'''i®22  N#  '‘•l®12  -  “2®11>I 

which  is  the  same  as  Equation  (61)  on  Page  42. 


This  equation  is  in  the  form 


^  +  E{<l>)  •  u  -I-  T{(l))  =  0 


(108) 


which  has  the  general  solution 

u  =  |c  -  yre-^^'^dc^ 


where 


cot  (j) 


f^d<l)  =  -Jcot</)d0 


N  ,  N 

IT  ^^1^22  “  ^2^12^  "  IT  ^^1^12  ■  ^2^11^ 


Then, 


After  integrating  the  simple  functions. 


Hence, 


u  =  jyi[N 


#’“'i®22  ■  ''2^:2’  ‘^e*“'i®i2  ■  ''2 


w  =  -ucot# 


^11^3  — T  C>sin^ 

1 1  sin  0  ( 


(109) 


(109a) 


(109b) 


(110) 


For  a  uniformly  symmetrical  distributed  loading,  p 
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d</) 

sin  <j) 


sin  ^ 


Similarly, 


(111) 


w  = 


-u  cot  (j) 


+  S 


12 


Specializing  for  spherical  shells  (Figure  14) 


=  N 


(112) 

NR 

w  =  -ucot  (ji  +  —  (S^^  +  S^^) 

where,  the  moduli  of  compliance  S  ,  S._  and  S„_  are  functions  of  <j) . 

11  1  ^  ^  ^ 

This  is  the  consequence  of  the  variation  of  the  filament  pattern  along 
the  meridian. 

At  the  juncture  with  the  insert,  the  boundary  of  the  shell  requires 

u  =  0  @  4)  =  (j)^ 

Therefore, 


(112a) 
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Edge  Effect  Deformation  of  a  Shell  of  Revolution  for  an  Upper  Edge 


In  what  follows,  the  asymptotic  solution  (Ref.  2)  is  used 


'r  1  r  * 

=  -—fx  d? 

=  /[c  cos?  +  C  sin?]e  ^  d? 

fO  f  X. 


=  — [C  (sin?  -  cos  I )  -  C-(sin?  +  cos?)]  +  C„ 

2^/1  ^  ^ 


at  ?  =  ? 

0 


*  m 

^ 


By  the  same  argument  on  the  order  of  magnitudes  of  the  terms. 


-w. 


m 


^  m  ^  m 

^0^20  ,  ^  ^  __!o_ 

r^  3  sin  4> 


which  yields  =  0 


Hence, 


*  1  f.^m 

*  '  -Ehl”e 


■'“XL,  ■-20 


r„_e  (sin?  +  cos?) 


*  "20L^m  ..,,.m 

w  = 


EhK-‘'”rL,  l=in(?  -r„)  +  cos(E  -r„)|e 

5  -  Sq 


-(?-?q) 


Similarly, 


(113) 


2r, 


1/4 


Eh  ^hr^Q  \  e 


-(r-?„)  . 

e  0  sin(£ 


?  =  r 


?  ) 
^  o' 
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(94) 


?  =  ?  - 


[3(1  -  i/^)] 


1/4 

(<^  -  </>q) 


Total  deformations  will  then  be 

*  .  m 
w  =  w  +  w 


m 

u  =  u 


X 


X  + 


m 

X 
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APPENDIX  --  STRESS  TRANSFORMATIONS 


Transformations  Relations  for  Anisotropic  Stresses 

Under  a  rotation  of  axes  from  x,  y,  z  to  P,  the  stress  components 
transform  in  accordance  with  the  following  relations 

O’  =  I-  I  ■  A  - 1 

p.q  ip  jq  ij 

where 

P,  q  =  €.  n,  ? 

i,  j  =  X,  y,  z 

I  =  direction  cosines 
'^Pg 

For  any  coordinate  axes; 


X 

y 

Z 

c 

4x 

4 

r] 

a 

px 

1 

ny 

L 

?y 

Expanding  the  expression  (A-1)  for  p  *  |  and  q  *  C 
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A -2 


XX  y?  yy  z§  zz 

+  2^  /  a  +2/  /  O'  +2/  I 

^■^xC  yC  xy  ^xg-^z?  xz  '^yC^zg  yz 


Similarly 


.2  ^  ^2  ^  2 

cr  -  i  Q  +  S  a  +  I  a 

m  XT)  XX  yr]  yy  ^zri  zz 


+  2i  io-  +21  la  +21  la 
xr}  yrj  xy  xrj  zr7  xz  yrj  zr;  yz 


o'j.j.  =  +  l^  +  l^  ^(y 

??  X?  XX  y?  yy  z?  zz 


+  2/  la  +21  la  +21  la 
^x^^y?  xy  ^xrz?  yz  z?  yz 


For  p  *  €;  q  *  r7 


/  /  O'  +i  L  o  +  I  - 1  a 

x§  xP  XX  yC  yp  yy  z§  zp  zz 

(^C^yp  ^  ^yc4n)°'xy  ^  ^ 


+  (/  /  +  I  /  W 

\V  zn  z?V/  yz 


For  p  *  C;  q  *  ? 


I  f.i  ^<y  +  /  »/  ^o  +  I  A  „o 

x§  X?  XX  yg  yl  yy  z§  z?  zz 


■^(4§4?^^y€4?j°'xy  ^  (^x?4?  ^  4§'4?)'^xz 
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For  p  *  n;  q  ®  ? 


7  ~  I  L  ^  /  o’  ■(■/  /  o’ 

nt  XX  ^ynj^  yy  zr/  Z?  zz 


A-7 


+  i  i  J 

la 

yh  x^l 

1  xy 

+  /  /  V,' 

1 

zn  y?, 

/  yz 

x)7  zr]  xH  xz 


)" 


Stress  Transformations  for  an  Orthotropic  Sheet 


Analysis  of  an  orthotropic  sheet  reduces  to  a  two  dimensional  problem 
because  the  stress  and  strain  components  normal  to  the  sheet  surface,  are 
assumed  to  vanish. 

This  is  equivalent  to  rotation  about  the  z-axis  (Figure  15).  Therefore, 


X 

y 

Z 

COS  (J 

cos  (f  -  u) 

=  sin  to 

0 

n 

cos  (f  +  w) 
=  -  sin  u 

cos  U) 

0 

0 

0 

1 

prj 


2  2 

o  cos  w  "t  o  sin  w  +  2o  cos  w  sin  <»; 
XX  yy  xy 

2  2 

a  sin  0)  +  a  cos  w  -  2a  sin  tocos  w 
XX  yy  xy 


A-8 


-  a  cos  (osin  to 

XX 


+  •  ^  /  2  .  2  \ 
+  a  sin  tocos  to  +  cos  to  -  sm  to  ) ' 

yy  \  / 


xy 
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